Thermal entanglement in superconducting quantum-interference-device qubits 

coupled to cavity field* 
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In this paper, we investigate thermal entanglement in a superconducting-quantum-interference- 
device qubit coupled to a cavity field. We show that the entanglement can be manipulated by varying 
temperature and an effective controlling parameter B which depends on the external field and 
characteristic parameters of the system. We find that there exists a critical value of the controlling 
parameter Be. Under a fixed temperature, increasing B can increase entanglement in the regime of 
B < Be, while the entanglement decreases with increasing B in the regime oi B > Be- 
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For decades, quantum entanglement has been the focus 
of much work in the foundations of quantum mechanics. 
Beyond this fundamental aspect, creating and manipu- 
lating of entangled states are essential for quantum in- 
formation applications. Hence, quantum entanglement 
has been viewed as an essential resource for quan- 
tum information processing. In recent years, much at- 
tention has been paid to the presence of entanglement in 
condensed-matter systems at finite temperatures. The 
state of a typical condensed-matter system at thermal 
equilibrium (temperature T) is described by the den- 
sity operator p = e~^^ /Z where H is the Hamiltonian, 
Z = tre~^^ is the partition function, and (3 — 1/kT 
where k is Boltzmanns constant. The entanglement as- 
sociated with the thermal state p is referred to as thermal 
entanglement 0|- 

Superconducting Josephson devices (SJD's) are sim- 
ple but realistic and extensively studied solid systems. 
In fact, they may be thought of as two-state sys- 
tems realizing the elementary unit of quantum informa- 
tion, known as a q ubit. Moreover, Josephson devices 
[1 H nil II II iM mill 111 can be scaled up to a large 
number of qubits and their dynamics may be controlled 
by externally applied voltages and magnetic fluxes. Su- 
perconducting devices such as Cooper pair boxes, Joseph- 
son junctions, or superconducting quantum interference 
devices (SQUID's) have been thus proposed and used as 
basic elements for the practical realization of quantum 
gates. 

It is an interesting topic to study thermal entangle- 
ment in systems involving SJD's, such as SQUID's. As 
we know, so far no study on thermal entanglement of 
Josephson qubit systems has been reported. The pur- 
pose of this paper is to investigate thermal entanglement 
of such a Josephson qubit system which consists of a flux 
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qubit (an rf-SQUID) coupled to a single-mode quantized 
electromagnetic field of a resonant cavity. 

The system which we consider is an rf-SQUID coupled 
to a single-mode electromagnetic field of a high-Q res- 
onant cavity. The rf-SQUID is a superconducting loop 
interrupted by a Josephson junction with the junction 
capacitance C and Josephson coupling energy Ej. The 
Hamiltonian of the rf-SQUID Q is given by 



Hs ^^-Ej cos{2n^) 
26 00 
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where is an externally applied dc flux, (j)Q = h/2e 
the flux quantum, L the self-inductance of the loop, the 
charge Q on the junction capacitance and the flux </> in 
the loop are canonically conjugate operators obeying the 
commutation relation [0, Q] = ih. 

If the self-inductance is large, such that Ej > 
(/<Q/(47r^L) and the externally applied flux <j)x is close to 
<j>o/2, the last two terms in ^ form a double- well poten- 
tial near (jj = (/'o/2. Actually, the form of the potential in 
Hamiltonian can be tuned by changing the external dc 
magnetic flux (f>x applied to the loop to the case where the 
two lowest-energy wells are degenerate. This case occurs 
for (px = 4>{)I2. These two degenerate minima correspond 
to the clockwise and counterclockwise sense of rotation 
of the supercurrent in the loop. The two relative ground 
states are localized flux states, hereafter denoted \L) and 
\R). At low temperatures only the lowest states in the 
two wells contribute. Hence the rf-SQUID reduces to an 
effective two-state system 10] with the Hamiltonian 



Hs 



(2) 



where the diagonal term Tie = 2/c-\/6(/3 — l)(</'x — 1^0/2) 
is the bias, i.e., the asymmetry of the double well, Tit 
can be tuned by the applied flux (px- The off-diagonal 
term A describes the inter- well tunnelling amplitude that 
depends on the Josephson coupling energy Ej which may 
be manipulated by a dc flux (j)c in a dc- SQUID replacing 
the single Josephson junction in the rf-SQUID . 
The eigenvectors of Hamiltonian (|2J) arc given by 



h) - Af_[|i?) + (e+ V^^+A^)/A|L)], (3) 
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FIG. 1: Entanglement as a function of the effective controlling 
parameter for different values of temperature: T = 0.1 (star 
line), T = 0.2 (dashed line), and T = 0.3 (solid line). Here B 
is in unit of 10~^^ J s. 



1+) = N+[\R) + ie-Ve^ + A^)/A\L)], (4) 

where the normalization constants N± = (1 + [(e ± 
Ve^ + A2)/A])~-'-/^, and the corresponding eigenvalues 
are E± = {±h/2)Ve^ + . 

The Hamiltonian of the single-mode cavity field can be 
written as He = hwc{a^a + l/2) where ujc is the frequency 
of the cavity field. One may model the cavity field as an 
LC resonator with the resonant frequency tOc — 1/ y/LcCc 
where the capacitive parameter Cc depends on the cavity 
field frequency via the quantization volume and on the 
rf-SQUID geometry [l^. 

The flux qubit and the iC-resonator modelling the 
single-mode cavity field can be thought of as coupled to- 
gether inductively. The interaction Hamiltonian |l3 can 
be described by 

Hj = B{a + a'f)(-ecr^ + Acr^), 



(5) 



where the controlling parameter B depends on the cou- 
pling strength and the system characteristic parameters, 
it can be explicitly expressed as 



B 



k 



L\ 2lJcCc + A2 ' 



(6) 



where k is the fiux linkage factor with the order of 0.01. 
Hence, the total Hamiltonian can be written as 



H = Hs + Hc + Hi. 



(7) 



If we focus our attention on the initial condition \na) — 
\n)\a) such that no more than one matter-radiation ex- 
citation n + (jia + 1/2) with rig. = ±1 is present. Under 
this initial condition, it can be shown that most of the 
evolution occurs within the low-lying energy subspace of 
the free Hamiltonian Hq = Hg -\- He spanned by the four 
states Incr), with n = 0, 1 and a = —,+. In such a four- 
state subspace, the dynamics system is governed by the 



reduced Hamiltonian 



H 




(8) 
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FIG. 2: Entanglement as a function of temperature for differ- 
ent values of the controlling parameter: B = 2.0 (solid line), 
B = 2.5 (dashed line), and B = 3.0 (star line). Here B is in 
unit of 10"^^ J s. 

We now investigate the thermal entanglement in the 
combined system of the rf-SQUID and the cavity field. 
For the simplicity, we consider the symmetric case of e = 
0. In this case, the eigenvectors of Hamiltonian (jSJ are 
given by 



=(|o- 



|l-»,l^2) 



1 



(|0+) + |l-)), (9) 



m = iV3(|l+)-6+|0-)),|V4>-Af4(|l+>-6-|l+», 

where the normalization constants A^3 = (1 + fe^)^^/^, 
and Ni = {1 + fe2^)-i/2 ^ith b± = ih± \/WTf?)/B, 
and corresponding eigenvalues are given by 

El = {h-B)uJc, E2^{h + B)ujc 
Es = {h- Vb^ + h^)u;c, Ei = {h+ \/b^ 



(10) 



We use the entanglement of formation as the entan- 
glement measure for a state described by a density op- 
erator p. To calculate this measure, starting from the 
density matrix , one can introduce the product matrix R 
associated with the density matrix and its time-reversed 
matrix R = p(crf o'|)p*(cri cr|) Then concur- 

rence is defined as C = TOaa;{Ai — A2 — A3 — A4, 0} where 
the Xi are the square roots of the eigenvalues of R, in 
decreasing order. In this method the standard basis, 
|0+), |0— ), |1-|-), |1— ) must be used. The entanglement 
of formation is a strictly increasing function of concur- 
rence; thus there is a one-to-one correspondence. The 
amount of entanglement is given by entropy defined as 



E = Hil 



1 + Vl - C2 



(11) 
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where h(x) = —xlog2X — (1 — x)log2{l ~ x). 

For the composite system of the rf-SQUID and the 
cavity field in thermal equilibrium at temperature T, 
making use of the eigenvectors and eigenvalues given in 
expressions ^ and (|10|l . we can get the corresponding 
density operator p(T) = J2i=ii<^i\'^i){^i\)/ where = 
exp {—PEi){i = 1,2,3,4) and Z = J^t^i^i ^^'^ P^^'" 
tition function of the composite system with (3 — 1/ kT 
with k being the Boltzmann's constant. In the standard 
basis, {| + 0), I + 1), I — 0), I — 1)}, the density operator 
can be written as 



PiT) = ^ 



f u V 

X w 

w y 

\v u. 



(12) 



temperature, entanglement exhibits a peak, which cor- 
responds to the critical value of the controlling param- 
eter. Entanglement increases with increasing B in the 
controlling regime of i? < Be, while entanglement de- 
creases with increasing B in the regime oi B > Be] (ii) 
entanglement at the critical point of the controlling pa- 
rameter Be increases while the critical point may move to 
the left with increasing temperature; (iii) in the limit of 
_B — > oo, entanglement approaches zero which means that 
the composite system disentangled. This large B behav- 
ior can be obtained analytically. Indeed, when B ^ h, 
we have b± = ±1, Ei « E3 and E2 ~ In this 

case, ^/r2 is the largest one among the square roots of 
the eigenvalues of the product matrix R, and the con- 
currence C = max {—62/(61 -f 62), 0} — which implies 
that entanglement vanishes E = 0. 



where the nonzero matrix elements are given by 

U = -(62 + 61), = -(62-61), 



5^63 6^64 



l + b'l 1 + 
5+63 6-64 



y = 



63 



64 



1 + 61 l + ^'l 



l + b'l 1 + bl' 



(13) 



where we have set b± = (^h± \^ B^ + /B. 
Then the product matrix R becomes 



R 



f +v'^ 
w'^ + xy 2wx 



V 2m?; 



2wy w + xy 



u'^ + v^ 




,(14) 



which has the following eigenvalues 

n - {w + ^yz\ r2^{w-^)yz\ 

r3 = {u + vf/Z^, r^ = {u-vf/Z^. (15) 

In order to calculate the amount of entanglement, we 
have to know which one among the square roots of the 
eigenvalues of R is the largest. This depends on the val- 
ues of the controlling parameter and temperature. Hence, 
it is difficult to analytically to calculate the amount of 
entanglement, and we turn to numerical investigation. 
In Figure 1 we have plotted how entanglement between 
the rf-SQUID qubit and the cavity field varies with the 
controlling parameter B for different values of tempera- 
ture T . From Fig. 1 we can see that (i) for every given 



Finally we look to the dependence of entanglement on 
temperature. In Figure 2 we have plotted the variation of 
entanglement with temperature for different fixed values 
of the controlling parameter B in the regime below the 
critical value B^- Figure 2 indicates that for every fixed 
value of B entanglement decreases slowly with increasing 
temperature in the low temperature regime of T < 0.15 
while entanglement drops rapidly in the middle temper- 
ature regime of 0.15 < T < 0.35, and disentanglement 
appears in the high temperature regime of T > 0.35. 
From Fig. 2 we can also see that entanglement increases 
with enhancing the value of the controlling parameter B. 



In summary, we have investigated thermal entangle- 
ment in a SQUID qubit coupled to a cavity field. We 
have shown that the entanglement can be manipulated by 
varying temperature and an effective controlling parame- 
ter B which depends on the external field and character- 
istic parameters of the system. In general, thermal en- 
tanglement increases with increasing the controlling pa- 
rameter and decreases with increasing temperature. We 
have found that there exists a critical value of the control- 
ling parameter Be- Thermal entanglement between the 
SQUID and the cavity field exhibits different properties 
below and above the critical value. Under a fixed tem- 
perature, increasing B can increase entanglement below 
the critical value, while the entanglement decreases with 
increasing B above the critical value. We hope that the 
finding of the existence of the critical value of the control- 
ling parameter is useful to realize quantum information 
processing by using SQUID qubits. 
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